Chiral Kosterlitz-Thouless transition in the frustrated Heisenberg antiferromagnet 

on a pyrochlore slab 
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Ordering of the geometrically frustrated two-dimensional Heisenberg antiferromagnet on a py- 
rochlore slab is studied by Monte Carlo simulations. In contrast to the kagome Heisenberg anti- 
ferromagnet, the model exhibits locally non-coplanar spin structures at low temperatures, bearing 
nontrivial chiral degrees of freedom. Under certain conditions, the model exhibits a novel Kosterlitz- 
Thouless-type transition at a finite temperature associated with these chiral degrees of freedom. 



Magnetic ordering of geometrically frustrated antifcr- 
romagnets (AFs) has attracted interest of researchers in 
magnetism In geometrically frustrated AFs, spins 

usually sit on lattices made up of triangles or tetrahedra 
as elementary units, and interact antiferromagnetically 
with their neighboring spins. Intrinsic inability to simul- 
taneously satisfy all antiferromagnetic nearest-neighbor 
interactions on a triangle or on a tetrahedron necessarily 
leads to macroscopic frustration. This makes the spin 
ordering on these lattices a highly nontrivial issue. 

Recently, interest has been focused on the properties 
of the lattices consisting of corner-sharing triangles or 
tetrahedra, e.g., the 2D kagome lattice and the 3D py- 
rochlore lattice Q-Q> Due to the looser coupling among 
the frustrating units, these systems often remain para- 
magnetic down to low temperatures. Experimentally, 
however, many of the geometrically frustrated magnets 
regarded as typical kagome or pyrochlore AFs exhibit a 
phase transition at a low but finite temperature, quite 
often a spin-glass (SG)-like freezing transition 

One of the best studied geometrically frustrated AFs is 
the S = 3/2 Heisenberg kagome AF, SrCrGaO (SCGO) 
HI). Experimentally, this material exhibits a SG-like 
transition at a finite temperature T = Tt as in many 
other geometrically frustrated AFs, although Tf is con- 
siderably lower than the Curie- Weiss temperature of this 
material due to the strong geometrical frustration 
In spite of extensive experimental and theoretical efforts, 
the true nature of this SG-like transition of SCGO has 
remained elusive. Although SCGO has been regarded 
for some time as a typical model compound of the 2D 
kagome AF, Monte Carlo (MC) simulations performed 
for the antiferromagnetic Heisenberg model on the 2D 
kagome lattice failed to reproduce the SG-like transition 
as experimentally observed in SCGO J|,[l(J, suggesting 
that the modeling of SCGO as the kagome AF might be 
inadequate in capturing some essential aspects of this ma- 
terial. In fact, the underlying lattice structure of SCGO 
is not of a pure (single-layer) kagome lattice, but rather, 
of a kagome sandwich, or a "pyrochlore slab" 
The structure of the lattice is illustrated in Fig.l: It 
consists of two 2D kagome layers which sandwiches the 
sparse triangular layer in between. Note that this lattice 
is obtained by slicing the 3D pyrochlore lattice along the 
(111) direction into the slab geometry. Inelastic neutron- 



scattering measurements have indicated that the neigh- 
boring slabs are magnetically well separated along the 
c-axis 11 . Hence, in modeling SCGO, one may safely 
neglect the inter-slab interaction and consider the 2D 
Heisenberg model on a pyrochlore-slab lattice. 
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FIG. 1. A pyrochlore-slab lattice, (a) The unit cell of the 
lattice consists of two corner-sharing tetrahedra illustrated by 
the solid lines: (b) The lower kagome layer consisting of the 
sites 1, 2 and 3, denoted by the solid circle. The open circle 
denotes an apical site of tetrahedron (site 4) . The upward tri- 
angles, shaded in the figure, correspond to the bottom planes 
of tetrahedra, and are grouped into three types A, B and 
C, each forming the three inter-penetrating triangular sublat- 
tices. 
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The purpose of the present Letter is to study the order- 
ing properties of the antiferromagnetic classical Heisen- 
berg model on a pyrochlore slab by means of MC sim- 
ulations, and to examine whether some new features 
which are different from those of the well-studied pure 
kagome Heisenberg AF would arise, possibly due to the 
tetrahedron-based structure of this lattice. In particular, 
we pay attention to the possible "chiral" properties of the 
model. "Chirality" is a multi-spin quantity representing 
the sense or handedness of the local non-coplanar spin 
structures induced by spin frustration. It is defined for 
three neighboring Heisenberg spins as a pseudo-scalar, 
X = Si ■ S-2 x S3, so as to give a nonzero value if the three 
spins make non-coplanar configurations but vanish oth- 
erwise. In the case of the pure kagome Heisenberg AF, 
it has been known that the spin structure selected at low 
temperatures is a coplanar one with the vanishing chiral- 
ity ||[l0]]. In sharp contrast to this, we show in the case 
of the pyrochlore-slab Heisenberg AF that the spin struc- 
ture stabilized at low temperatures is a non-coplanar one 
sustaining the nontrivial chirality. These nontrivial chiral 
degrees of freedom exhibit a novel thermodynamic phase 
transition at a finite temperature without accompanying 
the order of Heisenberg spins. 

The unit cell of the pyrochlore-slab lattice may be 
taken as two corner-sharing tetrahedra containing seven 
sites, as illustrated by the solid lines in Fig. 1(a). These 
unit cells containing seven sites, each numbered from 1 
to 7, are arranged forming the 2D triangular lattice of 
spacing 2d, d being the lattice constant of the kagome 
layers. In Fig. 1(b), we show the lower kagome layer con- 
sisting of the sites 1, 2 and 3. Note that the upward 
triangle in Fig. 1(b) corresponds to the bottom plane of 
tetrahedron, while the downward triangle is not a part 
of any tetrahedron. These upward triangles are further 
grouped into three types, A, B and C, each forming tri- 
angular sublattices of spacing 

We consider the classical Heisenberg Hamiltonian, 

yj = jjjSj ■ Sj + j s m ■ s n , (1) 

(ij) {ran) 

where the antiferromagnetic interaction Jy on the 
kagome layers is assumed to work between the nearest- 
neighbor (nn) and the next-nearest-neighbor (nnn) pairs, 
with Jij = J\ and J2, respectively, while J is the an- 
tiferromagnetic nn interaction between the kagome lay- 
ers and the triangular layer. The variable Si is a three- 
component unit vector at the i-th site. While we have 
made simulations for several choices of J2/J1 and J'/Ji, 
we report here on our results for the case of antiferro- 
magnetic nnn coupling J 2 > 0, because our model turns 
out to exhibit most interesting transition behavior in this 
case. 

Since the present model is highly frustrated showing 
very slow relaxation at low temperatures, we combine the 
heat-bath MC method with the temperature-exchange 



technique to facilitate efficient thermalization [ |12f . Sim- 
ulations are made for a pyrochlore slab with N — 7 x N s 
(N s = LxL) spins with L =6, 12, 18, 24 and 30. Periodic 
boundary conditions are employed. 

Now, we present our Monte Carlo results for the par- 
ticular case of J2/J1 = 0.5 and J' = J\. In Figs. 2(a) and 
(b), we show the specific heat and the uniform suscepti- 
bility, respectively. The specific-heat data exhibit double 
peaks, at T = T pl ~ 0.29 Ji and at T = T p2 ~ 0.09 Ji. 
Their size dependence reveals that the peak heights even- 
tually saturate with L, suggesting that both peaks are 
non-divergent ones. Meanwhile, the susceptibility data 
exhibit a clear cusp-like anomaly at T j J\ ~ 0.085 close 
to the lower specific-heat peak. 
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FIG. 2. The temperature and size dependence of (a) the 
specific heat and (b) the uniform susceptibility, per spin. The 
insets are magnified views of the low-temperature region. The 
arrow indicates the transition point. 

One generally expects that, in a fully isotropic Heis- 
neberg model in 2D like our model, there should be no 
long-range order (LRO) of Heisenberg spins, nor a finite- 
temperature transition occurring in the spin sector. We 
have confirmed this expectation by calculating various 
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order parameters and Binder ratios associated with the 
spin order, e.g., the so-called q — and y/3 x V3 or- 
ders on the kagome layer |^,|l0|. We also find that the 
higher specific-heat peak is correlated with the develop- 
ment of the q — short-range order of Heiscnberg spins 
on the kagome layers. In the present model with the an- 
tiferromagnetic J2, the q = mode is selected over the 
Vox vo mode, in contrast to the pure kagome Heisenberg 
AF where the %/3 x %/3 mode is selected HQ . 
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summation is taken over all upward triangles and (■ • •) 
represents the thermal average. The data shown in the 
inset of Fig. 3 (a) indicate that in the T — > limit x tends 
to a nonzero value, meaning that the spins belonging to 
tetrahedron form locally non-coplaner structures sustain- 
ing nontrivial chirality. The appearance of the nontrivial 
chiral degrees of freedom in the present model is in sharp 
contrast to the behavior of the pure kagome Heisenberg 
AF. 

In order to detect the possible ordering of these chiral- 
ities, we calculate the ferrimagnetic (or staggered) chiral 
order parameter m x f, defined by 
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This quantity gives a nonzero value if the chirality ex- 
hibits a ferrimagnetic order with -\/3 x periodicity on 
three sublattices A, B and C in Fig. 1(b). The associated 
chiral Binder ratio may be defined by 
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FIG. 3. (a) Log-log plot of the L-dependence of the chi- 
ral susceptibility associated with the v3 X v3 ordering at 
several temperatures. The inset represents the temperature 
dependence of the mean local amplitude of the chirality for 
L = 24. The L-dependence of the data is negligible here, 
much smaller than the size of the symbols, (b) The tempera- 
ture dependence of the chiral Binder ratio for various L. The 
inset is a magnified view of the transition region. The arrow 
indicates the transition point. 

The chirality-related quantities are shown in Fig. 3. We 
define here the local chirality, \ii f° r the three spins 1, 
2 and 3 located at each upward triangle on the lower 
kagome layer. In order to detect the local non-coplanarity 
of the spin structures, we compute the mean local ampli- 
tude of the chirality by \ 2 = (1/^Vs) YliiXi )> where the 



In Fig. 3(b), we show the calculated chiral Binder ra- 
tio. For smaller L, g x f for various L tend to cross at 
T = T c ~ 0.082Ji, but for larger L they tend to merge 
at T < T Cl signaling the occurrence of a phase transi- 
tion of the chirality. The estimated transition temper- 
ature T c /Ji — 0.082(2) is in rough agreement with the 
susceptibility-cusp temperautre estimated above, and is 
slightly below the lower specific-heat peak temperature 
T p 2- In fact, there is no apprecialbe anomaly in the 
specific heat just at T = T c : See the inset of Fig. 2(a). 
A merging behavior of the Binder ratio, without the 
discernible specific-heat anomaly just at T c but with a 
nondivergent peak slightly above T c , suggests that the 
observed chirality transition might be of the Kosterlitz- 
Thouless (KT)-type. 

In Fig. 3(a), we show on a log-log plot the L-dependence 
of the chiral susceptibility N s {m x f) associated with the 
ferrimagnetic chiral order. While the data at higher tem- 
peratures exhibit the behaviors characteristic of the dis- 
ordered phase, bending down toward some finite values, 
those at T < T c lie on straight lines, exhibiting the be- 
havior expected for the KT-like phase with algebraically- 
decaying correlations. The estimated slope of the plots, 
which should be equal to 2 — r\ with r\ being the critical- 
point decay exponent, turns out to be 77 = 0.24 ± 0.01 at 
T = T c and gradually decreases with decreasing temper- 
ature below T c . This again indicates that the observed 
chiral transition at T — T c is the KT-type transition. 
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Since the chirality here is an Ising like quantity sit- 
ting on the triangualr lattice, there is a close similarity 
between the chirality ordering of the present model and 
the ordering of the 2D Ising model on the triangular lat- 
tice. Indeed, when the triangular Ising model possesses 
the antiferromagnetic nn interaction and ferromagnetic 
nnn interation, the model is known to exhibit a v3 x yo 
order via the KT-type transition characterized by the ex- 
ponent r) — 1/4 [13|-|l6|]. Since the chirality ordering of 
the present model with J 2 > is essntially the staggered 
one, it would be no surprise that the chirality ordering 
here is of the KT-type. To the authors' knowledge, our 
present finding is the first case of chiral KT transition 
without accompanying the spin order. 

It should be noticed that the triangular Ising AF with 
the ferromagnetic nnn interaction exhibits another phase 
transition with further decreasing temperature, into the 
low-temperature phase with a finite LRO. fp^|-^6| The 
exponent ri at this second transition is believed to be 
1/9 p^ , p^ , p^7t . In view of this, we search for this sec- 
ond transition into the long-range-ordered state in our 
present model, but with negative result. Presumably, 
severe frustration inherent to the present model might 
hinder the onset of the true LRO. 

Although we have so far considered the case where the 
inter-plane interaction J' is equal to the intra-plane nn 
interaction J\ , such equality is not expected in general in 
real experimental systems. In order to examine the possi- 
ble effect of varying the inter-plane coupling, we also per- 
formed similar calculations for the J' = 0.5 Ji case, to get 
qualitatively similar results as the J' = J\ case. We thus 
conclude that the </' value is irrevant to the chiral KT 
transition observed here. Further detailed results of our 
MC simulations, including the results for other choices of 
the nnn interaction J2, will be published elsewhere [ pj] . 
For the cases of the vanishing and ferromagnetic nnn in- 
teraction J2 < 0, we find that, although the local spin 
structure is non-coplanar sustaining nontrivial chirality, 
the chiral KT transition does not take place there. 

Finally, we wish to refer to the possible experimental 
implication of our present results. We emphasize that the 
present model is expected to capture essential geometri- 
cal ingredients of SCGO, and indeed, seems to accout for 
some of the experimental features, e.g., the occurrence of 
a finite-temperature phase transition above all, the cusp- 
like anomaly observed in the susceptibility , or the ex- 
istence of a broad specific-heat peak slightly above T c . 
However, we mention that some other experimental fea- 
tures remain unexplained. For example, experimentally, 
the transition is characterized by the negative divergence 



of the nonlinear susceptibility and the onset of magnetic 
irreversibility, e.g., a notable difference between the field- 
cooled and zero-field-cooled magnetizations Jjj . We could 
not reproduce these features in our simulations. 

Although further work is clearly required to fully 
account for the experimental observation on SCGO, 
our present finding of the chiral KT transition in the 
pyrochlorc-slab Heisenberg AF might hopefully give a 
new perspective in the studies of SCGO, and geomet- 
rically frustrated AFs in general. 

The numerical calculation was performed on the Hi- 
tachi SR8000 at the supercomputer center, ISSP, Uni- 
versity of Tokyo. 
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